Hole polaron formation and migration in olivine phosphate materials 
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By combining first principles calculations and experimental XPS measurements, we investigate 
the electronic structure of potential Li-ion battery cathode materials LiMP04 (M=Mn,Fe,Co,Ni) 
to uncover the underlying mechanisms that determine small hole polaron formation and migration. 
We show that small hole polaron formation depends on features in the electronic structure near the 
valence-band maximum and that, calculationally, these features depend on the methodology chosen 
for dealing with the correlated nature of the transition-metal d-derived states in these systems. 
Comparison with experiment reveals that a hybrid functional approach is superior to GGA+U in 
correctly reproducing the XPS spectra. Using this approach we find that LiNiP04 cannot support 
small hole polarons, but that the other three compounds can. The migration barrier is determined 
mainly by the strong or weak bonding nature of the states at the top of the valence band, resulting 
in a substantially higher barrier for LiMnP04 than for LiCoP04 or LiFeP04. 

PACS numbers: 



I. INTRODUCTION 

LiFeP04 is an olivine structured material with many 
properties that make it attractive for usage as a cath- 
ode in Li-ion rechargeable batteries. It is composed of 
inexpensive starting materials, and the existence of a 
stable end compound, FeP04, allows for full withdrawal 
of Li ionsii However, the intrinsic electronic conductiv- 
ity of (T = 1.8xl0~* S/cm is prohibitively \ow^ so that 
special processing techniques, such as nanostructuring, 
coating with carbon or doping with supervalent cations 
must be employed in order to make the material suit- 
able for electrochemical cycling. '^^^ The underlying mech- 
anism by which these techniques are successful is still 
hotly debated)2rJ^ but whatever it is, the conductivity 
can be reliably raised by a factor of 10^-10^, render- 
ing LiFeP04 an excellent, cycleable cathode material. 
LiMnP04, LiCoP04, and LiNiP04 are all isostructural 
to LiFeP04 and have higher voltages, both calculated 
and measurediii^"— , suggesting they could be used as 
cathode materials with even more overall energy than 
LiFeP04. Unfortunately, these materials have similar 
or even poorer intrinsic conductivities r^"— measured at 
cr(LiMnP04)< 10-1" S/cm, cr(LiCoP04)-- lO'^ S/cm, 
and (T(LiNiP04)'^ 10^^ S/cm. To compound the prob- 
lem, techniques to raise these numbers have been only 
moderately successfullSr— . The conductivities can be 
raised by 10^-10^, but problems of strong capacity fade, 
perhaps associated with the still high resisitivity, con- 
tinue to hamper their usefulness as practical cathode ma- 
terials. 

The olivine phosphates are now widely understood to 
be wide band-gap materials that exhibit polaronic rather 
than band-like transport.— Their electrochemically 
active center is the transition-metal ion with formal va- 
lency 2"*" in the stoichiometric compound, that becomes 



S"*" upon withdrawal of a Li ion and associated electron 
or when, e.g., Li vacancies are created during synthcsis^^'. 
The resulting localized hole and consequent contraction 
of the surrounding O ions are together known as a small 
hole polaron. In order to move through the crystal, the 
small polaron must hop from one transition-metal site to 
another. The barrier to this hopping creates the acti- 
vated transport seen in experiment j^i^i^^ 

In this work, we use first principles density-functional 
theory (DFT) and x-ray photoemission spectroscopy 
(XPS) to examine the effect of the electronic structure 
of LiMP04 (M=Mn,Fe,Co,Ni) on small hole polaron for- 
mation and migration with the goal of understanding 
the differences in measured conductivities and poten- 
tially pinpointing techniques to raise them. We find that, 
calculationally, the formation of small hole polarons de- 
pends on which approximation to the exchange correla- 
tion functional is used. By comparing finely detailed fea- 
tures of the density of states (DOS) in our computational 
results to careful XPS measurements, we can determine 
which methodology provides the best agreement with the 
true electronic structure. We determine that the use of 
a hybrid functional that allows localization of electrons 
on both transition-metal and oxygen sites is necessary 
to reproduce the observed binding energies. Using this 
functional, we calculate formation energies and migra- 
tion barriers for small hole polarons in each of the com- 
pounds. We find, in good agreement with experiment, 
that LiFeP04 hole polarons have the highest tendency 
to hop, followed by LiCoP04 and LiMnP04. We find 
that small hole polaron formation is highly unlikely in 
LiNiP04, which may also become electronically unstable 
with any significant delithiation. The migration barriers 
depend strongly on the lattice distortion which, in turn, is 
determined by details of strong and weak bonding states 
formed with surrounding oxygen ions. 
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II. METHODOLOGY 
A. Computational 

Our calculations employ the projector augmented wave 
mcthodr^i21 as implemented in the VASP code.^^'^^ For 
the exchange correlation potential, we use the Perdew- 
Burke-Ernzerhof version of the generalized gradient ap- 
proximation (PBE-GGA)^ the GGA+U correctional 
(U values were taken from the calculations of Zhou 
et alM-) , and the Heyd-Scuseria-Ernzerhof hybrid func- 
tional (HSE06)<^2'2i The structures are fully relaxed 
within each type of calculation. The nudged elastic band 
(NEB)'^'^ method was used to calculate the migration bar- 
rier of a polaron moving from one site to another. 

To obtain our calculated binding energy spectra, we 
projected out the character of each ion from the DOS, 
and adjusted its intensity according to the Scofield cross- 
sections for x-rays^ at 1487 eV, relative to the intensity 
of the transition-metal ion. The adjusted partial DOS 
are then summed and convolved with a Gaussian func- 
tion to approximate the temperature broadening seen in 
experiment. This spectra is compared to the experimen- 
tal binding energies measured by XPS. 

B. Experimental 

LiMnP04 was prepared by a hydrothermal method us- 
ing urea hydrolysis as the source of hydroxide ions.'^'' 
MnS04-H20, Li2S04-H20, 85% H3PO4 and urea were 
dissolved in enough water to make a 20 mL clear so- 
lution. The molar ratio, Li:Mn:P04:urea, was 1:1:1:1.5. 
The solution was heated to 200° C in sealed. Teflon-lined, 
45 mL autoclave for 15 hours. The product was isolated 
by centrifuge and dried under vacuum at 80° C. Phase 
purity was determined via x-ray powder diffraction us- 
ing a Rigaku Ultima III x-ray diffractometer in a parallel 
beam geometry. Lattice constants were determined via 
Rietveld refinement of the parallel beam x-ray data us- 
ing the Riqas program. Lattice constants of the obtained 
LiMnP04, in spacegroup Pnma, were a = 10.456 A, b — 
4.7503 1, c = 6.1006 A, unit cell volume = 302.99 A^. 
LiFeP04 was prepared similarly using FeS04-7H20 in 
place of MnS04-H20. The lattice constants of the ob- 
tained LiFeP04 were a =10.3180 1, 6=4.6901 A, c = 
5.9960 1, unit cell volume = 290.17 A^. LiCoP04 sam- 
ples were prepared via a citrate complexation route<^ 
Co(OH)2, LiH2P04, and citric acid, 1, 1.01, 1.02, molar 
ratio, respectively, were mixed into deionized water un- 
til all solids were dissolved. The resulting solution was 
evaporated to dryness via a microwave oven. The result- 
ing dried mass was removed, ground lightly with mortar 
and pestle and heated in air at a rate of 10° C min"^ to 
600° C and the reactant mixture was held at this temper- 
ature for 12 h. Lattice constants were a = 10.1950 A, b 
= 5.9179 A, c = 4.6972 A for a unit cell volume of 283.40 
A^. LiNiP04 was obtained from the solid state reaction 



of LiH2P04 with Ni(0H)2. The starting materials were 
mixed in a stoichiometric ratio via a mortar and pestle. 
The mixture was heated at 325° C for 12 hours, 500° C 
for 8 hours, 600° C for 12 hours and 700° C for 12 hours 
with intermittent grinding. All heating steps were done 
in air. Lattice constants were a = 10.0319 A, b — 5.8522 
A, c = 4.6779 A, unit cell volume of 274.63 A^. 

XPS spectra were acquired using a Kratos Axis 165 
x-ray photoelectron spectrometer operating in hybrid 
mode, using monochromatic Al ka radiation (1486.6 eV) 
at 220 W. The powder samples were attached to the sam- 
ple holder using double sided conductive copper tape, 
charge neutralization was required to minimize sample 
charging. The spectrometer was at 5 x 10^^ Torr or 
lower throughout data collection. Survey spectra and 
Valence band high resolution spectra were collected with 
pass energies of 160 eV and 20 eV respectively. All data 
were calibrated to the hydrocarbon contamination peak 
at 284.8 eV. 



III. RESULTS 

A. Electronic structure and hole polaron formation 

Since the creation of a hole polaron involves removing 
an electron from the valence-band maximum ( VBM) , the 
nature of the electronic states at and near the VBM is 
crucial in determining if the polaron can be formedi^ We 
have identified the following two criteria that are neces- 
sary to computationally establish small hole polarons in 
the olivine phosphate materials: (1) There must be a 
narrow band at the VBM that produces a sharp peak 
in the DOS, with a finite gap separating it from the 
broad continuum of states below it; and (2) the com- 
position of this band must be predominantly transition- 
metal-derived, with 0-derived character comprising less 
than approximately 1/3 of the total. Furthermore, some 
corrective method must be applied to deal with the cor- 
related transition-metal-derived d states. In uncorrected 
GGA, these are overpenalized for localization and imme- 
diately delocalize throughout the crystal, whether or not 
the two aforementioned criteria are satisfied. In Fig. [1] 
we plot the total DOS, as well as the partial DOS (shaded 
region) that comes from transition metal character (re- 
maining character may be assumed to be oxygen) for the 
four olivine phosphate compounds. As can be seen in 
Fig. [T]and Table HI both criteria are filled for LiFeP04, 
regardless of which exchange correlation potential is used 
(GGA-I-U or HSE06). For LiMnP04 and LiCoP04, only 
the first criterion is filled, while the second depends on 
which methodology is applied. The difference between 
the two results is significant: if oxygen character domi- 
nates at the VBM, the hole created upon electron with- 
drawal will be delocalized, whereas if transition-metal 
character dominates, small hole polarons can form and 
there will be a typical redox center. Furthermore, the 
predominance of oxygen character at the VBM means 
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that withdrawal of an electron produces oxygen ions that 
differ significantly from their preferred 2~ valency, elec- 
tronically destabilizing the compound and resulting in a 
tendency towards oxygen recombination and release as 
gas42ii^ For LiNiP04 neither correction produces the 
necessary features. In the GGA+U calculation, the Ni 
character has been pushed well below the oxygen states, 
leaving almost none at the VBM. Although some Ni re- 
mains at the VBM in HSE06, the dominant contribution 
is still oxygen. Therefore, small hole polaron formation 
cannot be achieved. This indicates that electronic trans- 
port in this compound must proceed via another mech- 
anism, and that electronic instability during charging is 
likely. 
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FIG. 1: DOS plots for LiMP04 (M=Mn,Fe,Co,Ni) in the an- 
tiferromagnetic state. The zero of the energy is set to the 
highest occupied state. "Up" spin states are plotted on the 
positive y-axis and "down" spin states on the negative y-axis. 
Shading corresponds to projected transition-metal character 
with majority "up" (dark) or majority "down" (hght) in the 
antiferromagnetic ordering pattern. The left hand side of the 
Figure shows the GGA-I-U generated DOS, while the right 
side shows the HSE06 generated DOS. Considerable differ- 
ences in oxygen admixture in the valence band are noticeable, 
except for M=Fe. 



B. XPS binding energies 

To establish which of our methodologies, GGA+U or 
HSE06, better reproduces the actual electronic structure 
in the olivine phosphates, and therefore to determine 



TABLE I: Percentage of transition-metal character in band 
nearest to the VBM in LiMP04 (M=Mn,Fe,Co, Ni) for three 
different exchange correlation schemes. The character not 
attributable to the transition-metal at the VBM comes almost 
exclusively from oxygen. 





Transition metal character at VBM 
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GGA+U 


HSE06 


LiMnP04 


71% 


39% 


60% 


LiFeP04 


88% 


85% 


89% 


LiCoP04 


87% 


8% 


78% 


LiNiP04 


94% 


9% 


34% 



whether or not small hole polaron formation is possible, 
we compare both to measured XPS spectra in Fig. [2] For 
all four compounds, the peak positions are almost per- 
fectly captured by the IISE06 calculation, whereas the 
GGA+U spectra are significantly compressed resulting 
in binding energies that are too low compared to experi- 
ment. Additionally, there are obvious differences near the 
VBM between the two methodologies and in all cases, 
the experimental spectra are much better matched by 
HSE06 than by GGA+U. The spectra as a whole are 
better represented at both high and low binding energies 
using IISE06 in comparison to GGA+U. 

Galculated peak heights are somewhat distorted com- 
pared to measurement for all four compounds and in both 
methodologies. This can be attributed to the fact that 
we have to project out the separate atomic characters 
from the DOS to apply the Scofield cross-section param- 
eters, and not all states are equally localized within the 
radii used for projection. This is most problematic for 
oxygen character which tends to be the most delocalized 
and therefore most missed in the projection. Oxygen 
has a small experimental photoionization cross-section so 
that missed character in the DOS remains uncorrected 
and produces exaggerated peak heights compared to ex- 
periment. The two peaks at highest binding energy, at- 
tributable to strongly bound P and O states, are stronger 
than in experiment in all four compounds. Since P has a 
high cross-section and O a low one, small errors in the ra- 
tio of the two will cause the peaks to be artificially high. 
The same effect is operative in the LiNiP04 spectrum 
where Ni and O are very strongly hybridized. The delo- 
calized oxygen character leads to the exaggerated peak 
at -3 eV. 

Because an electron is withdrawn from the very top 
of the valence band, the character of the states there is 
extremely important for hole polaron formation. If tran- 
sition metal character dominates, charge localization is 
possible and hole polarons will form. If oxygen charac- 
ter dominates, charge will delocalize and polarons will 
not form. The very good match between IISE06 binding 
energies and those measured in experimental XPS, indi- 
cates that the character of these states has been well- 
captured by this methodology, as strong shifts in the 
character would also show up as shifts in peak positions 
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FIG. 2: A comparison of experimental XPS measurements of the binding energies to DFT calculations with two different 
exchange-correlation approximations. In all cases the overall width of the spectrum and peak alignments are far better repro- 
duced using the HSE06 methodology. 



as can be seen in the GGA-I-U calculations. All spectra 
have two sharp peaks at high binding energy attributable 
to heavily bonded P and O states (the two leftmost ar- 
rows in Fig. (2). The next peak, located between -4 eV 
and -2 eV represents not a single kind of state, but a 
broad continumm of states composed of transition- metal 
and oxygen bonding and antibonding states which give 
rise to moderate secondary structure in the spectra, es- 
pecially in in the calculations. A peak near the valence 
band maximum has previously been identified as the lo- 
calized state necessary for small hole polaron formation 
in LiFeP044i In all of our spectra, we see a similar low 
energy peak peak in all of our spectra, but a comparison 
with the unbroadened and character-resolved DOS (Fig 
[T]) shows that in LiNiP04, it is not transition-metal dom- 
inated, nor is it even separated from the continuum of 
states below. In all three other compounds, this peak rep- 
resents a separate, energy-localized and transition-metal 
dominated state in which polaron formation is possible. 



Temperature broadening mainly obscures the separation 
between this "polaron state" and the continuum, but 
close inspection of the experimental spectrum does allow 
the interpretation that LiNiP04 has only three clearly 
distinguishable peaks (plus a small shoulder) due to a 
merged first and second peak, whereas the other com- 
pounds have four. Comparison with calculation makes 
this interpretation the most reasonable. The lowest en- 
ergy peak in LiNiP04 is more closely associated with 
the second lowest peak in the other compounds and is 
not suitable as a "polaron state". We conclude that the 
VBM of LiNiP04 is dominated by oxygen character, that 
no gap between the initial peak and the rest of spectrum 
exists and thus, small hole polarons will not form during 
delithiation of LiNiP04 and oxygen will likely be evolved. 
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FIG. 3: Energy of polaron migration as calculated at steps 
along a linearly interpolated path between two calculated po- 
larons. Comparison with a previous calculations is provided 
for reference. 



C. Small hole polaron conduction 



To understand how small hole polarons move through 
the ohvine structure, we used a method similar to that 
of Refs. .42,4J} to calculate the migration barrier for a po- 
laron hopping from one transition metal site to another. 
First we remove an electron from a given transition-metal 
site in a sixteen formula unit supercell of LiMP04. We 
fully relax all the ionic positions to achieve the oxygen 
distortion around the hole, along with any other nearby 
distortions that occur. In a separate calculation, we re- 
move an electron from an adjacent transition-metal site 
and again perform a relaxation. We find that the distor- 
tion of the lattice is negligible beyond a single unit cell, 
and is mainly constrained to the first shell of neighbors 
surrounding the hole, confirming the small polaron des- 
ignation. We use a linear interpolation of the two struc- 
tures and calculate the energies along a series of inter- 
mediate positions, using the NEB method. The highest 
energy we consider to be the polaron migration barrier. 
The positive (hole) charge density in this scenario "fol- 
lows" the distortion, hopping wholly or partially from 
one transition-metal site to the next without ever oc- 
cupying an intermediate space. The migration barrier, 
therefore reflects mainly the energy required to distort 
the lattice. As can be seen in Fig. [31 the resulting num- 
bers are 0.33 eV for LiMnP04, 0.20 eV for LiFeP04, 
and 0.23 eV for LiCoP04, following the trend of mea- 
sured conductivities in these materials. The numbers for 
LiMnP04 and LiFeP04 are shifted slightly higher, but 
otherwise in good agreement with Ref. |43| using the same 
exchange correlation potential. 



IV. DISCUSSION 



A. Effect of exchange-correlation corrections 



In a Hartree-Fock calculation, the subtraction of the 
Hartree and exchange terms produces an exact cance- 
lation of the electron self-interaction energy. When an 
approximation to the exchange term is made in DFT, 
this exact cancellation is compromised and leftover self- 
interaction contributes to the eigenvalues and total en- 
ergy of the system.''^ This is especially problematic for 
systems with localized bands. Both the GGA-I-U and 
HSE06 methods are useful to counteract the erroneous 
self-interaction term, but they are different in their ef- 
fects on the electronic structure. The GGA-I-U method 
applies a correction to the correlated bands (transition- 
metal d-bands for the purposes of this article) only. Since 
all non-correlated bands are unaffected, one effect is that 
the occupied (down-shifted) d bands are brought closer 
in energy to the oxygen bands below them, compared to 
GGA alone4^ The hybridization between the O p bands 
and metal d bands depends inversely on the energy dis- 
tance between them, which systematically decreases with 
the increasing value of the transition-metal ion. Conse- 
quently, the smaller energy separation in GGA-I-U causes 
a large increase in the amount of oxygen mixed into 
the transition-metal-derived bands at the valence band 
maximum (See Table U and Fig. [T]). For LiCoP04 and 
LiNiP04, where the d states are already quite low in 
energy, the downshift is so dramatic that the d-derived 
bands are located within or even below the oxygen band 
complex resulting in a predominantly oxygen VBM char- 
acter, which does not accurately represent the experi- 
mental data. 

The HSE06 method, on the other hand, incorpo- 
rates exact exchange along with the GGA approxima- 
tion to the exchange-correlation potential, with each be- 
ing used in a certain region of real space. '^'^ The exact 
exchange portion does not require any discrimination 
between "correlated" and "non-correlated" orbitals and 
therefore shifts both metal-derived and 0-derived states 
downward, though not equally. Compared to GGA-fU, 
the mixture of O character into the highest energy va- 
lence states is greatly reduced. Comparison with the 
XPS spectra show that this reduced hybridization bet- 
ter reproduces the experimental binding energies. For all 
but LiNiP04, this results in a sharp peak with mainly 
transition-metal character at the VBM and a gap to the 
continuum of states below it. For the Ni-based com- 
pound, oxygen character at the VBM is reduced from the 
GGA-I-U value, but is still too high to allow small hole 
polaron formation. The good match between experiment 
and calculation shows that this is due to hybridization 
resulting from the legitimately low energy position of the 
Ni d bands, and not due to computational artifacts. 
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TABLE II: Comparison of M-oxygen bond lengths for 
LiMP04 (M=Mn,Fe,Co) at M sites with a locahzed hole po- 
laron (M"^"*") and without a hole polaron (M^''"), obtained from 
supercell calculations. Contraction around the hole occurs in 
all cases, although a Jahn- Teller distortion expands two bonds 
in LiMnP04). 



LiMnP04 


M-0 bond lengths 


Ave 


Mn^+ 
Mn^+ 


2.13 2.13 2.15 2.25 2.27 2.27 
1.93 1.95 2.00 2.00 2.34 2.34 


2.20 
2.09 


LiFeP04 






Fe^+ 
Fe'+ 


2.06 2.06 2.12 2.21 2.25 2.25 
2.01 2.01 2.02 2.07 2.13 2.14 


2.16 
2.06 


LiCoP04 






Co^+ 
Co3+ 


2.04 2.05 2.05 2.11 2.16 2.18 
1.94 1.94 2.00 2.10 2.13 2.14 


2.10 
2.04 



B. Lattice changes during polaron formation 



crystal. Such distortions are also likely to be detrimental 
to the structure during repeated electrochemical cycling, 
especially at high rates. 
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LiFePa 




The lattice contraction around the removed electron 
depends strongly on the occupation of the metal d states 
prior to withdrawal. In the olivine structure, the transi- 
tion metal sits in a quasi-octahedral MOg environment. 
The crystal field and ligand interactions split the 5 d 
states into a lower three-fold degenerate manifold (t2g) 
and an upper doubly degenerate manifold (eg). The 
states are strongly exchange split such that the ground 
state is high spin in all cases. For LiMnP04, the hole 
state is in a eg state, while for LiFeP04 and LiCoP04, 
the hole sits in a t2g state (See Fig. [H). The anti-bonding 
Cg states are heavily mixed with oxygen - the result of 
strong cr-type bonding between metal and oxygen states, 
whereas the weak- or non-bonding t2g states are much 
more lightly mixed with oxygen. Removing an electron 
from the anti-bonding states increases the bond strength 
and the surrounding octahedron contracts strongly. Re- 
moving an electron from the non-bonding states results 
in a much weaker shift of oxygen toward the transition- 
metal ion. This can be seen in Table HI] where the average 
Mn-0 bond distance changes much more than Fe-0 or 
Co-0. Lattice distortion around the Mn'^+ ion is compli- 
cated by the fact that creation of a hole leaves a single 
electron in the doubly degenerate eg complex, stimulat- 
ing a Jahn- Teller distortion which partially counteracts 
the overall contraction of the MnOg octahedra. The two- 
long/four-short pattern typical of octahedral Jahn- Teller 
distortions is combined with an overall shrinking of all 
Mn-0 bonds to produce four dramatically shorter bonds 
and two slightly longer ones. Thus, the average contrac- 
tion in bond length (5%) belies the actual magnitude 
of the distortion: a 9% contraction of four bonds and a 
3% expansion of two others. The average contraction of 
bond lengths in LiFeP04 is 4.3% and in LiCoP04 is 2.7%. 
The strong distortion in the Mn-based material will dis- 
rupt small hole polaron migration (electronic conduction) 
since large local changes must propagate through the 



FIG. 4: Charge density of the small hole polaron in phosphate 
olivine compounds. Most of the positive (hole) charge resides 
on the transition-metal, but some charge is also on surround- 
ing oxygens, particularly in the case of LiMnP04 for which 
the hole sits in a strong bonding, heavily oxygen-hybridized 
state. Below each polaron density plot is a schematic show- 
ing the energy levels and fillings in the (approximately) oc- 
tahedral environment of each compound. Strong-bonding eg 
and weak-bonding t2g states are show with the majority spin 
states (left) widely exchange split from the shaded minority 
spin states (right). The heavy, hollow arrow indicates the 
position of the removed electron or hole polaron. 



V. SUMMARY 

Our investigation of four olivine phosphate compounds 
reveals that small hole polaron formation depends on the 
existence of a narrow, isolated band at the VBM that 
contains predominantly transition-metal character. The 
relevant electronic structure is sensitive to the particular 
approximation used for the exchange-correlation poten- 
tial. GGA+U pushes the transition-metal bands down 
much further relative to the oxygen bands than does 
HSE06, thereby increasing the hybridization and mix- 
ing more oxygen into the valence band. A comparison 
with experimental XPS spectra shows that only HSE06 
properly captures the true electronic structure for these 
materials at the VBM. Using this methodology, we find 
that LiMnP04, LiFeP04, and LiCoP04 can all support 
small hole polarons, while LiNiP04 cannot due to a pre- 
dominantly oxygen-derived valence band that delocalizes 
the hole upon electron removal. The hole polaron state in 
LiMnP04 is within the strong-bonding eg complex and 
additionally supports a Jahn- Teller distortion, resulting 
in an overall stronger local deformation of the lattice and 
therefore higher migration barrier than in LiFeP04 or 
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LiCoP04 where the hole polaron occupies a weak- or 
non-bonding state. The order of barriers is LiFeP04 < 
LiCoP04 < LiMnP04, in agreement with experimental 
measurements of conductivity. 
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